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Abstract 

A fully automated iterative design method has been 
developed by which an airfoil with a substantial amount 
of natural laminar flow can be designed, while 
maintaining other aerodynamic and geometric 
constraints. Drag reductions have been realized using 
the design method over a range of Mach numbers, 
Reynolds numbers and airfoil thicknesses. The thrusts 
of the method are its ability to calculate a target N- 
Factor distribution that forces the flow to undergo 
transition at the desired location; the target-pressure-N- 
Factor relationship that is used to reduce the N-Factors 
in order to prolong transition; and its ability to design 
airfoils to meet lift, pitching moment, thickness and 
leading-edge radius constraints while also being able to 
meet the natural laminar flow constraint. The method 
uses several existing CFD codes and can design a new 
airfoil in only a few days using a Silicon Graphics IRIS 
workstation. 

Nomenclature 

A relaxation factor for N-Factor design method 

C p pressure coefficient 

c ; section lift coefficient 

c pitching moment coefficient 


* Graduate Research Scholar Assistant, AIAA Student 
Member 

f Professor of Engineering and Applied Science, Asso- 
ciate Fellow of AIAA 

$ Research Engineer, Senior Member of AIAA 
§ Research Engineer 

Copyright © 1996 by the American Institute of Aeronautics and 
Astronautics, Inc. No copyright is asserted in the United States under 
Title 17, U.S. Code. The U.S. Government has a royalty-free license 
to exercise all rights under the copyright claimed herein for govern- 
ment purposes. All other rights are reserved by the copyright owner. 


N N-Factor 

s chordwise distance from the leading-edge 

x non-dimensional distance along chord of airfoil 

Greek 

X relaxation factor used to meet lift constraint 

5 relaxation factor used to meet pitching moment 

constraint 

Subscript 
a analysis 

cp control point 

des design 

j airfoil station 

stag stagnation point 

T target 

to I tolerance 

u upper surface 

Superscript 

i iteration index 

Introduction 

For many years, reducing drag has been a major 
issue in airplane design. Reductions in drag allow 
airplanes to operate more efficiently by using less fuel, 
which results in reduced operating costs and smaller, 
quieter engines. The design of airplanes with bigger 
payloads and longer ranges is also possible when drag is 
reduced. 

There have been many attempts to reduce airplane 
drag. Since the 1930’s, there has been a great interest in 
designing airfoils and wings for natural laminar flow 
(NLF) to reduce viscous drag. In addition, supercritical 
wings were developed to reduce the wave drag on 
airplanes. The use of winglets on airplanes today also 
demonstrates an attempt to reduce airplane drag. This 
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paper presents a design method for reducing drag 
through the design of NLF airfoils. 

This new design method 1 combines proven CFD 
analysis and design codes with a new target pressure 
design technique. A flowchart of the method is shown 
in figure 1. Starting from an initial airfoil, the flow 
solver, which is an Euler solver - ’ ’ coupled with a 
turbulent boundary layer method 5 , is used to calculate 
the pressure distribution of the starting airfoil. Then, a 
laminar boundary layer solver 6 is used to calculate the 
boundary layer velocity and temperature profiles. These 
profiles are then used by COSAL 7 , the stability analysis 
code, to calculate the N-Factors of the current airfoil. 
Using the current pressures and N-Factors, the Target 
Pressure Design module then calculates a target pressure 
distribution that increases the amount of laminar flow 
over the airfoil. 

Once the target pressures are known, the CDISC 
airfoil design method 8 then iteratively designs a new 
airfoil which possesses this target pressure distribution. 
The target pressures for the design are then modified in 
the Modify Target Pressures to Enforce Constraints 
module while using the airfoil design method such that 
the aerodynamic constraints are achieved 1 . After 
designing a new airfoil, the boundary layer profiles are 
then calculated once again by the laminar boundary 
layer solver so that the stability analysis code can 
calculate the N-Factors of the new airfoil. If all of the 
design constraints are met, then the method terminates. 
Otherwise, the method will continue to iterate until the 
prescribed number of iterations is achieved. 

The Target Pressure Distribution 

After choosing an initial airfoil, the pressure 
distribution of the airfoil is calculated using the Euler 
solver, and the N-Factor envelope of the airfoil is 
calculated using the laminar boundary layer solver and 
the stability analysis code. Then, from the current 
pressure and N-Factor distributions, a target pressure 
distribution is calculated that meets all of the 
aerodynamic constraints. The module labeled Target 
Pressure Design in figure 1 is responsible for calculating 
this new target pressure distribution. Figure 2 shows a 
detailed flowchart of how the Target Pressure Design 
module calculates the new target pressures. Each of the 
components in this flowchart will now be discussed. 

The Target N-Factor Distribution 

Once the analysis N-Factor distribution has been 
calculated by the stability analysis code, a target N- 
Factor distribution must be prescribed, as demonstrated 
by the first module in figure 2. This target N-Factor 
distribution must force the flow to undergo transition at 
the desired location. 



Figure 1 . A flowchart of the NFF airfoil design method 

In order to calculate the target N-Factor 
distribution, four control points (x cp j , x cp 2 , x cp 3 and 
x cp 4 ) are specified. The first control point is located at 
the point where the analysis N-Factors first become 
greater than an N-Factor level, IV l . Ahead of the first 
control point, the analysis N-Factors are kept as the 
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Figure 2. A flowchart of the Target Pressure Design 
module 


target N-Factors. Beyond this point, the other control 
points and their N-Factor levels are used to calculate the 
target N-Factors. This is done by connecting each 
consecutive control point with a straight line, then 
smoothing the curve using a polynomial fit. 

Figure 3 shows a typical target N-Factor 
distribution and the control points used to calculate 
these target N-Factors. In this figure, x cp 0 = 0.58, 
x , = 0.60 , x . = 0.65 , N , = 3 , N , = 8 , 

cp, 3 ’ cp, 4 ’ c/7, 1 ’ cp, 2 ’ 

N , = 10, and A , = 15 . Ahead of x ~ , the target 
N-Factor distribution forms a buffer region above which 
the N-Factors aren’t allowed to grow so that the flow 
will not undergo transition at slightly off-design 
conditions. Then, between x , and x , , the N- 
Factors are allowed to grow rapidly to force transition at 
60% chord for this case. 



Figure 3. A typical target N-Factor distribution with 
control points 

Moreover, the third control point is placed at the 
location of desired transition, and a nominal transition 
N-Factor of 10 is used. Actually, the transition N-Factor 
could have been as low as eight or as high as 15. For 
this purpose, the N-Factors were allowed to grow 
rapidly between the second and fourth control points. 
This reduces the uncertainty in the transition location by 
increasing the N-Factor gradient in the transition region. 

Extrapolation of Analysis N-Factors 

The stability analysis code can calculate N-Factors 
only when the boundary layer is attached. This causes 
problems if the boundary layer of the current airfoil 
separates ahead of the fourth control point, x cp 4 . The 
second module in figure 2 is included to remedy this 
situation. 
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Suppose the boundary layer of the current airfoil 
separates at airfoil station k, say at 30% chord. In 
addition, suppose that the fourth control point is located 
at airfoil station / , say at 65% chord. This indicates that 
analysis N-Factors exist only to 30% chord, while the 
target N-Factors have been specified back to 65% chord. 
In the N-Factor design method, which will be discussed 
in the next section, both analysis and target N-Factors 
must be known ahead of the fourth control point. 


For this reason, a method was developed to 
extrapolate the analysis N-Factors from station k to 
station 1 . The method extrapolates the N-Factors based 
on the pressures between the two stations. At any 
station j between k and Z, the analysis N-Factors are 
calculated as 


N, 


N 


N. 


N j - 1 + C 


7-1 7-2 


_ r 

p,a,j-\,u p,a,j-2,u 


x (C 


c . , ) 

p, a,j - 1 ,u' 


(i) 


Where C p,a,j,u 

pressure coefficient at station j . 


P, a,j, it 

represents the upper surface analysis 


The N-Factor Design Method 

After calculating the target N-Factor distribution 
and extrapolating the analysis N-Factors, the N-Factor 
design method calculates a target pressure distribution 
that moves the analysis N-Factors toward the target N- 
Factors. This method, which is shown as the third mod- 
ule in the flowchart in figure 2, is based on a target pres- 
sure-N-Factor relationship that was developed. In this 
method, the change in pressure coefficient required at 
airfoil station j to move the analysis N-Factor at j , N - , 
toward the target N-Factor at j , N r ■ , is 
AC . 

— Sd = AAA. (2) 

Ax. J 

J 

where 


AC . = C T . -C (3) 

p,] p, T.j, u p, a,], u ' ’ 

A *i = \ X j~ X j-l\ (4) 

AM. = N T .- Nj (5) 

In these equations, A is a relaxation parameter 
(typically 0.50), C p T ■ u is the upper surface target 
pressure at station j , C p a . u is the upper surface 
analysis pressure at station j, and Ax. is the grid 
spacing. 


In order to maintain a smooth and continuous target 
pressure distribution, once A C p . has been calculated at 
station j , this change in C p is applied to all of the sta- 
tions downstream of j as well. In doing this, the flow 
downstream of j has been changed and, as a result, it is 
necessary to correct the analysis N-Factors downstream. 
This is done by adding A M. to each of the analysis N- 


Factors aft of j . Equation 2 can then be rewritten to 
account for the corrections made downstream. The tar- 
get pressure coefficient that results at an airfoil station j 
is 

p.T.j.u v p.a.j.u p,j-\,u ’ 

+ AAx j (N Tj - (N j + AN j _ l )) (6) 

With j = 1 corresponding to the stagnation point, this 
relation is valid from j = 2 to j = l , where / is the air- 
foil station representing the fourth control point, x 4 . 
The boundary conditions at the stagnation point are 

C T . = C , = C , (7) 

p,l,l,u p,a,\,u p, stag x 7 

N t j = Mj = 0 (8) 

This method is used only to design the upper 
surface target pressures since the lower surface target 
pressures are only modified to meet the geometric 
constraints. 


The Pressures in the Recovery Region 

Once the target pressures have been calculated 
ahead of station l , the target pressures in the recovery 
region are calculated, as demonstrated by the fourth 
module on the flowchart in figure 2. These target 
pressures are determined based on the analysis pressures 
in the recovery region of the initial airfoil that was 
analyzed and a linear pressure coefficient distribution 
that is eventually used to modify the pressure 
distribution to meet the pitching moment constraint. 

In order to calculate the target pressures in this 

region, two intermediate pressure distributions are 

determined. The first intermediate pressure distribution, 
C p j | , is determined by linearly scaling the pressures of 
the initial airfoil, C p . 0 . (The linear scaling method 
used is described in Appendix A of reference 1.) As a 
result, the first intermediate pressure distribution is cal- 
culated using the relation 
C . , = C „ + 

p,J.\ p.n. 0 

c - c 

. P.T.l.u p,n,0 / n n \ , ni 

r -r y p.j.0 ^ p.n, o’ 

L p.1,0 L p.n,0 

where n is the trailing-edge station. This expression is 
valid from j = l to j = n . 


The second intermediate pressure distribution, 
C . , , is determined by adding a linear pressure coeffi- 
cient distribution to C . A . As a result, 

p.j.0 


C 


P.j, 2 


C pJ.Q + o.30 


x. 

-^C 


p, 0.7 c 


( 10 ) 


where C p Qlc is the magnitude of the linear pressure 
coefficient distribution at 70% chord. This equation is 
also valid from j = 1 to j = n . 
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Finally, the target pressures in the recovery region 
are calculated by taking a weighted average of the first 
and second intermediate pressure recoveries. (The 
weighted averaging technique is discussed in Appendix 
B of reference 1.) The upper surface target pressures in 
the recovery region then become 


C. 


p, T,j. u 


! lr 

-X P’j - 1 


— -C ■ , 

_v, P.J, 2 


(ID 


Equation 1 1 is used from j = I to j 


Meeting the Aerodynamic Constraints 

The upper surface target pressures that have been 
calculated must now be modified to meet the lift and 
pitching moment constraints. This is shown in the fifth 
module on the flowchart shown in figure 2. 


In order to modify the target pressure distribution to 
meet the aerodynamic constraints, the upper surface 
target pressure distribution is divided into three distinct 
regions. The leading-edge region extends from the 
leading-edge of the airfoil to the last station where the 
N-Factors are zero, which is represented by station m . 
The center region extends from station m + 1 to / , the 
location of the fourth control point. This is the region 
where the N-Factors are growing most rapidly. Finally, 
the recovery region extends from station l+l to n , the 
trading-edge. 


A method was developed to modify the target pres- 
sures in such a way that the upper surface NLF would 
not be disturbed while meeting the aerodynamic con- 
straints. A detailed flowchart of this process is shown in 
figure 4. Once c, and c m have been calculated, the lift 
constraint is checked. If |c, , - cj > c, , , , then AC 

| /, des l\ l, tol ’ p, c 

is calculated and the target pressures in the three regions 
are modified. With 


AC 


2 Xs (c 


/, des 


■C,) 


p, c 


( 12 ) 


the new target pressures in the center and leading-edge 
regions are, respectively, 

J+ i 


C T . — C T ■ + AC 

p.T.j.u P.T.J.U p, c 


(13) 


*->, T,j, u C p T i u 


(14) 


Cp, T,j, u ~ Cp t, 1, u + 

c' + 1 - c' , 

p, 1, m, u p, 1 , 1, u 

c' T - c' , 

p, l,m, u p, 1 , 1 , u 

where i denotes the old target pressures, i + 1 denotes 
the new target pressures, X is a relaxation factor (typi- 
cally 0.50) and 


= II- 

k = 2 


(15) 



Figure 4. A flowchart of how the upper surface target 
pressures are modified to meet the aerodynamic con- 
straints 
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Once the target pressures in the leading-edge and center 
regions have been calculated, the method of the 
previous section is used to calculate the new target 
pressures in the recovery region. As shown in the 
flowchart in figure 4, this process is iterated until 

|0, des 0| — 0, tol ' 

Once the lift coefficient meets the desired tolerance 
c, fol (typically 0.01), a new value of C p 07c is calcu- 
lated to modify the target pressures to meet the pitching 
moment constraint. The new value of C ,, n is calcu- 

p, 0.7 c 

lated as 


4 + 1 


24 c; (c 


■H _ r . . m,des~ C J 

P - a7e P • Q Jc ( 1 -x t ) (8.r, + 1) 


(16) 


where i denotes the old value of C and c is a 

/7,0.7 c ^ 

relaxation factor (typically 0.50). Once the new value 
of C p 07c has been calculated, the method of the 
previous section is used once again to calculate the 
target pressures in the recovery region. The target 
pressures are then modified again to meet the lift 
constraint. This method is iterated until 


u l,des O 


— C,| <c 


/, tol 


and 


, des 


<C 


tol ’ 


where 


c j is the pitching moment tolerance (typically 0.01). 
The development of this entire method can be found in 
reference 1. 


Once the upper surface target pressures have been 
modified to meet the aerodynamic constraints, the upper 
and lower surface target pressures are modified to meet 
the leading-edge radius constraints. Then, after 
modifying the upper surface target pressures to once 
again meet the aerodynamic constraints, the lower 
surface target pressures are calculated. These lower 
surface target pressures are determined by linearly 
scaling the upper surface target pressures. Once this is 
done, the upper surface target pressures are once again 
modified to meet the aerodynamic constraints. The 
processes of modifying the target pressures to meet the 
leading-edge radius constraint and calculating the lower 
surface target pressures will not be discussed here. 
They have been discussed in detail in reference 1 . 

Results 

The NLF airfoil design method that has been 
described in the previous sections will now be used to 
design several airfoils for a variety of flow conditions 
and constraints. 


Airfoil for a General Aviation Application 

The first airfoil that is presented was designed for 
the same flow conditions and constraints for which the 
NLF(1)-0414F airfoil 9 - 10 was designed. The NLF(l)- 
0414F has already been tested in a wind-tunnel. This 
airfoil was designed to have 70% chord NLF on both 
surfaces at a Mach number of 0.40, a Reynolds number 


of 10 million and a lift coefficient of 0.40. In order to 
have NLF to 70% chord, the N-Factor at 70% chord was 
to be less than nine on both surfaces. In addition, the 
airfoil was to be 14% thick and have a trailing-edge 
cruise flap with the pressures in the recovery region 
being determined using an energy method. When the 
boundary layer had the most energy, the pressure 
gradient was to be the most adverse. Then, as the 
boundary layer loses energy, the pressure gradient was 
to become less adverse so that turbulent separation 
would be avoided at the design condition. 

The NLF(1)-0414F was designed iteratively using 
the NYU code 11 - 12 - 13 to calculate the pressures of each 
intermediate airfoil and the SALLY code 14 to perform 
the stability analysis. Modifications were then made by 
hand to the airfoil geometry to give the airfoil the 
desired characteristics. 

Starting with the NACA 64J-212 airfoil, a new 
airfoil was designed by the NLF airfoil design method 
that would meet the same constraints as the NLF(l)- 
0414F. Since the NLF airfoil design method does not 
have the capability to use an energy method to prescribe 
the pressures in the recovery region, the pitching 
moment coefficient was constrained by the NLF airfoil 
design method to be -0.08 on the new airfoil. In 
addition, the airfoil thickness at 85%, 90% and 95% 
chord were constrained by the NLF airfoil design 
method to prevent the trailing-edge from being too thin. 

Figure 5 shows a comparison of the new airfoil and 
the NLF(1)-0414F. The airfoils designed by the two 
different techniques appear to be very similar. Since the 
leading-edge radius was not constrained in either 
design, there is a slight difference between the two 
airfoils in the leading-edge region which is reflected in 
the mid-chord region. If the leading-edge radius had 
been constrained, perhaps the airfoils would be even 
more similar. 

Figure 6 shows a comparison of the pressure 
distributions of the two airfoils. In general, these 
pressure distributions also appear to be very similar. 
There is a little difference in the leading-edge pressures 
due to the difference in the leading-edge radius of the 
two airfoils. The pressures in the mid-chord region 
compare very well although they are slightly off-set. 
This may be due to the fact that the NLF airfoil design 
method was able to maintain the thickness of the new 
airfoil to be exactly 14% chord, while the maximum 
thickness of the NLF(1)-0414F is actually 14.3% chord. 
The bumps in the pressure distribution on the lower 
surface of the new airfoil in the recovery region are due 
to the thickness constraints imposed at 85%, 90% and 
95% chord. 
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Figure 5. A comparison of the NLF(1)-0414F and the 
new airfoil 



X 


Figure 6. A comparison of the pressures of the NLF(l)- 
0414F and the new airfoil at M = 0.40, Re = 10 million, 
and C| = 0.40 

Figure 7 shows a comparison of the upper and 
lower surface N-Factor distributions of the two airfoils. 
The N-Factors at 70% chord on the new airfoil are less 
than nine, as was prescribed. The target N-Factor 
distribution discussed previously (see figure 3) was not 
used in this design since that type of distribution was not 
prescribed in the design of the NLF(1)-0414F. In figure 
7, the N-Factor distributions of the new airfoil possess 
the same general shape as the N-Factor distributions of 
the NLF(1)-0414F. In order to meet the lower surface 
N-Factor constraint, the lower surface of the airfoil was 
designed by inverting the airfoil and designing it at 
-c, . and -c . . This was done since the N-Factor 

l, des m, des 

design method discussed previously can only prescribe 
the target pressures on the upper surface of the airfoil so 
that the N-Factor constraint is achieved. 



X 


Figure 7. A comparison of the N-Factor distributions of 
the NLF(1)-0414F and the new airfoil at M = 0.40, Re = 
10 million, and C| = 0.40 

The design of the NLF(1)-0414F airfoil took 
approximately one year, while the design of the new 
airfoil using the NLF airfoil design method took only 
two days on a Silicon Graphics Indigo2 workstation 
with an R4000 processor. 

Airfoil for a Glider Application 

As the second application of the NLF airfoil design 
method, the glider airfoil presented in reference 1 will 
be redesigned to improve the laminar flow 
characteristics of the lower surface. The glider airfoil 
was designed for 65% chord upper surface NLF at a 
Mach number of 0.10, a Reynolds number of 3 million, 
a lift coefficient of 0.30, and a pitching-moment 
coefficient of -0.06. In addition, the airfoil was to be 
15% thick with a leading-edge radius of 1.4% chord. 
Front and rear spar constraints were also enforced. 

In redesigning this airfoil, 65% chord NLF is 
desired on both surfaces while maintaining the lift and 
pitching moment constraints, and the maximum 
thickness and leading-edge radius constraints. It was 
assumed that the front and rear spar thicknesses could 
not be constrained while maintaining the desired 
amount of NLF on the lower surface. As a result, the 
front and rear spar constraints were removed for the 
current design. 

Figure 8 shows a comparison of the glider airfoil 
from reference 1 and the redesigned airfoil. The 
redesigned airfoil meets the maximum thickness 
constraint of 15% chord and has a leading-edge radius 
of 1.38% chord. Notice that the redesigned airfoil is 
thicker than the glider airfoil in the region from 50% to 
80% chord. This increase in thickness was the 
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Figure 8. A comparison of the glider airfoil from refer- 
ence 1 and the redesigned airfoil 

mechanism by which a more favorable pressure gradient 
was imposed on the lower surface so that laminar flow 
could be increased. This can be seen in figure 9, which 
compares the pressure distributions of the two airfoils. 



X 


Figure 9. A comparison of the pressure distributions of 
the glider airfoil from reference 1 and the redesigned 
airfoil at M = 0.10, Re = 3 million, and C| = 0.30 

Figure 10 shows a comparison between the upper 
and lower surface N-Factor distributions of the 
redesigned airfoil and the target N-Factor distributions. 
In this plot, a nominal N-Factor of 13.5 was used to 
estimate the location of transition. On both surfaces, the 
N-Factor distributions meet the target N-Factor 
distributions. As a result, NLF was achieved on both 
surfaces to 65% chord. The plot also shows the lower 
surface N-Factor distribution for the glider airfoil of 
reference 1. The upper surface N-Factor distribution of 
the glider airfoil is not shown in this plot since it did not 
change much during the design. 



X 


Figure 10. The upper and lower surface N-Factor distri- 
butions of the redesigned airfoil at M = 0.10, Re = 3 mil- 
lion, and Cj = 0.30 

Airfoil for a Commuter Aircraft 

As the final example, the commuter airfoil in 
reference 1 was redesigned to increase the amount of 
NLF on the lower surface, since a laminar flow 
constraint was only implemented on the upper surface in 
the design of this airfoil. This airfoil was designed at a 
Mach number of 0.60, a Reynolds number of 20 million 
and a lift coefficient of 0.40. The airfoil was to be 12% 
thick, with a leading-edge radius of 1% chord. In the 
current design, it was desired to maintain 60% chord 
NLF on both surfaces with a pitching moment 
coefficient of -0.08. 

In the current design, repeated attempts were made 
to design a 12% thick airfoil with a 1% chord leading- 
edge radius, while meeting the NLF and aerodynamic 
constraints. In these attempts, the upper and lower NLF 
constraints were not being achieved while maintaining 
both the thickness and leading-edge radius constraints. 
As a result, the leading-edge radius constraint was 
released so that both of the NLF constraints could be 
achieved. Maintaining the leading-edge radius 
constraint would have required the maximum thickness 
of the airfoil to be larger than 12% chord. While 
maintaining the maximum thickness constraint, the 
leading-edge radius of the redesigned airfoil was 0.78% 
chord. All of the other constraints were achieved in this 
design. 

Figure 1 1 shows the redesigned airfoil and its 
pressure distribution at the design conditions. The 
favorable pressure gradients to 70% chord on both 
surfaces are responsible for the flow remaining 
laminar. Figure 12 shows the upper and lower surface 
N-Factor distributions of the airfoil at the design 
conditions. In this figure, an N-Factor of 10 was used to 
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estimate the location of transition. This critical value 
was varied from the previous example to show that the 
NLF airfoil design method works for various transition 
N-Factors. In this plot, the N-Factor distributions meet 
the target N-Factor distributions that were imposed in 
order to force the flow to undergo transition at 60% 
chord. 



Figure 1 1 . The redesigned commuter airfoil and its pres- 
sure distribution at M = 0.60, Re = 20 million, and q = 
0.40 

The upper surface target N-Factor distribution in 
this plot resembles the target N-Factor distribution that 
was discussed previously (see figure 3). The target N- 
Factor distribution in figure 12 forms a buffer region 
ahead of 50% chord where the N-Factors are not 
allowed to grow above eight. Then, from 50% to 70% 
chord, the N-Factors are allowed to grow rapidly to 
force the flow to transition at 60% chord. 



X 


Figure 12. The upper and lower surface N-Factor distri- 
butions of the redesigned commuter airfoil at M = 0.60, 
Re = 20 million, and q = 0.40 


Concluding Remarks 

An automated two-dimensional method has been 
developed for designing NLF airfoils, while maintaining 
several other aerodynamic and geometric constraints. 
The method has been shown to work for a range of 
Mach numbers, Reynolds numbers, and airfoil 
thicknesses. 

In order to develop this NLF airfoil design method, 
several existing CFD codes were coupled together. In 
addition, a process was developed for calculating a 
target N-Factor distribution that forces transition to 
occur at the desired location. Using this target N-Factor 
distribution, as well as the current analysis N-Factors 
and pressures, a method was also developed for 
calculating a target pressure distribution. Using this 
target pressure distribution, the current airfoil is 
redesigned to obtain a new airfoil that is closer to 
meeting the desired NLF, aerodynamic and geometric 
constraints. This method has been used to design a 
number of airfoils, with results shown for glider, general 
aviation and commuter applications. 

There are several possibilities for extension of this 
research. The method could be applied to bodies other 
than airfoils and wings, with possible applications 
including fuselages and nacelles. In addition, the 
method could be extended to the design of airfoils for 
supersonic applications. Since large sweep angles are 
needed for supersonic wings, crossflow instabilities 
would be a major issue. In these cases, boundary layer 
suction and blowing is often necessary to help reduce 
the crossflow disturbances. As a result, when extending 
the method to include supersonic designs, the method 
may also have to be modified to account for suction and 
blowing. 
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